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Sumaary 


The  characterisation  of  five  gas^liquid  chromatographic  stationary 
phases,  using  retention  data  obtained  by  Laffort  and  co-workers  for  240 
solutes,  has  been  carried  out  by  the  method  of  multiple  linear  regression 
analysis.  Relative  gas-liquid  partition  coefficients  for  as  many  solutes  as 
possible  were  correlated  against  various  combinations  of  parameters.  The 

best  general  equation  was  found  to  be  one  containing  the  exploratory  variables 

*  16  * 
u  ,6,a  ,0  ,  and  log  L  ,  Attempts  to  replace  “fi  by  the  dipole  moment 
2  m  ra  2 

2 

(u  )  were  not  very  success!  .1,  and  neither  were  attempts  to  use  refractive 

2 

index  functions  or  molar  refractions  in  combination  with  U  .  However, 

* 

replacement  of  dipole  moment  itself  was  more  successful,  and  led  to 

quite  good  equations  in  4,  0^,  and  log 

A  number  of  new  log  values  have  been  determined  experimentally,  using 

the  gas -chroma cographic  method . 
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Introduction 


1  2 

Over  the  past  few  years,  Abraham,  Doherty,  Kamlet,  Taft  and  co-workers  ’ 
have  constructed  equations  for  the  correlation  and  prediction  of  a  very  large 
number  of  physicochemical  and  biochemical  phenomena.  These  equations  are 
based  on  a  cavity  theory  of  solution,  in  which  the  process  of  dissolution  of  a 
solute  in  a  solvent  may  be  broken  down  into  a  number  of  hypothetical  steps: 

(i)  the  endoergic  formation  of  a  cavity  in  the  bulk  solvent,  (ii) 
rearrangement  of  solvent  molecules  round  the  cavity,  and  (iii)  the  exoergic 
interaction  of  the  solute  with  the  surrounding  solvent  molecules  after  the 
solute  has  been  inserted  into  the  cavity.  If  the  Gibbs  energy  change  in  step 

(ii)  is  zero,  or  very  nearly  zero  as  is  usually  assumed,  only  steps  (i)  and 

(iii)  need  to  be  modelled.  The  energy  of  formation  of  a  cavity  can  be  taken 

2 

as  proportional  to  the  solvent  cohesive  energy  density,  where  6^  is  the 

Hildebrand  solubility  parameter,  and  to  some  function  of  the  solute  size,  or 
volume,  ^2,'^  leading  to  a  term  units  of  energy.  Various 

solvent-solute  interactions  can  take  place  in  step  (iii).  If  both  solvent 
and  solute  are  dipolar,  a  term  in  arise,  where  is  the  solvent  or 

solute  di polarity.  Hydrogen- bond  interactions  will  also  be  set  up,  either 
between  a  solvent  acting  as  a  hydrogen-bond  base  and  a  solute  acting  as  a 
hydrogen- bond  acid,  or  between  a  solvent  acting  as  the  acid  and  the 


solute 

as  the  base , 

In 

these  two 

terms ,  and  3^  refer 

to  the 

solvent 

hydrogen- bond  acidity 

and 

basicity. 

and  and  32  to  the 

solute 

hydrogen-bond  acidity  and  basicity.  The  full  equation  for  the  correlation  of 
some  solubility  related  property,  SP,  is  then  given  by  the  multiple  linear 
regression  equation, 


We  denote,  as  usual,  solvent  properties  by  the  subscript  I  an;:  solute 
propf»rties  by  the  stjby,cript  2. 


(1) 


SP  =  A  +  Bn*u*  + 

Now  for  a  process  involving  a  series  of  solutes  in  a  given  solvent,  all 
the  solvent  parameters  in  equation  (I)  are  constant,  leading  to  equation 
(2).  For  solutes  that  are  aromatic  or  pol yha logenat ed ,  a  polarisabi 1 ity 

★ 

SP  =  C  +  s.n^  +  a,a^  +  +  M*V2  (2) 


correction  term  is  needed,  which  takes  the  form  ^2  ~  ^  aromatic  solutes, 
0.5  for  polyha logenated  solutes,  and  zero  for  all  other  solutes.  This  leads 
to  the  final  equation,  used  extensively  in  the  correlation  of  a  wide  variety 
of  phenomena  in  condensed  phases,  equation  (3).  An  early  application  of 


SP 


C  +  s.Ti^  +  a.a^  +  b. 


5^  +  ^*^2 


(3) 


equation  (3)  to  a  process  involving  a  gaseous  phase,  namely  the  solubility  of 
gases  and  vapours  in  polymers ,  revealed  a  possible  deficiency  in  that 
equation  (3)  contains  no  term  that  corresponds  to  solute-solvent  dispersion, 
or  van  der  Waals,  interaction.  An  alternative  equation  was  therefore  put 
forward ,  with  a  new  solute  parameter,  log  replacing  the  volume  terra 
V2.  This  new  parameter  was  defined  as  the  logarithm  of  the  solute  Ostwald 
solubility  coefficient,  L,  on  n-hexadecane  at  298K.^^ 

L  =  concentration  of  solute  in  solution 

concentration  of  solute  in  the  gas  phase 

Two  possible  equations  for  the  correlation  of  the  solubility  of  a  series 
of  gases  and  vapours  in  a  given  condensed  phase  are  therefore  equations  (3) 
and  (5),  and  we  set  out  to  investigate  the  use  cf  these  two  equations.  There 
are  several  ways  of  describing  the  solubility  of  gases  in  liquids,  but  in  view 


of  the  use  of  the  Ostwald  solubility  coefficient  in  equation  (5)  shall  define 


gas  solubility  through  equation  (4).  It  should  be  noted  that  h  is  actually 
the  same  as  the  gas-liquid  partition  coefficient,  K,  used  in  the  description 
of  gas-liquid  chromatography  (GLC),  The  coefficient  L  or  K  is  related  to  the 
specific  retention  volume  at  the  column  temperature,  V^,  through  equation  (6) 
where  o^  is  the  stationary  phase  density. 

SP  =  C  +  s.n*  +  ci.62  +  a.a^  +  +  l.log  l’'^  (5) 

L  (or  K)  =  (6) 


In  order  to  apply  equations  (3)  and  (5),  values  of  SP  (i.e.  log  L  or  Log 
Vq)  should  be  available  for  a  wide  selection  of  solutes  on  the  same  liquid 
phase.  By  far  the  most  convenient  and  accurate  method  of  obtaining  such  a 
series  of  SP  values  is  by  GLC  itself,  and  there  are  numerous  compi  lations  of 
such  data.  Probably  the  most  extensive,  carefully  measured,  values  are  those 
recorded  by  Laffort  and  co-workers^  for  240  solutes  on  5  stationary  phases, 
and  this  is  the  data  we  have  chosen  first  to  analyse.  Laffort  and  co-workers 
published  their  data  in  the  form  of  Kovat's  retention  indices,  defined  for 
isothermal  GLC  through  equation  (7). 


I(x)  =  100 


log  V^(x)-log  V^(Pn) 
1  og  ^  og 


+  lOOn 


(7) 


Here,  I(x)  is  the  retention  index  of  solute  x,  Vq  is  the  specific  retention 
volume,  and  (Pn)  and  ( Pn  +  1)  represent  n-alkanes  of  carbon  number  n  and  N  + 
1.  The  equation  used  to  calculate  their  retention  indices  is, 

log  t(x)  -  log 


I  =  100 


b 


+  1000 


(8) 


where  t(x)  and  corrected  retention  times  for  compound  x  and 

n-decane,  and  b  is  the  slope  for  n-aikanes.  From  equation  (6)  and  equation 
(8)  it  follows  that 

log  L(x)  =  b  +  log  L(P^q)  (9) 

and  hence  values  of  log  L(x)  may  be  calculated  for  the  240  solutes  relative  to 
log  L  for  n-decane.  These  relative  values  of  log  L(x)  -  log  LCP^q)  can  be 
used  in  equations  (3)  and  (5)  to  yield  exactly  the  same  coefficients  etc.  as 
would  be  found  with  the  absolute  values  of  log  L(x).‘'‘  The  constant  quantity 
log  will  be  subsumed  into  the  value  of  C. 

The  five  phases  studied  by  Laffort  and  co-workers  are  listed  in  Table  1, 
nd  the  240  solutes  are  given  in  Table  2.  In  Table  3  are  given  those  listed 
in  our  database,  i.e.  those  for  which  we  had  some,  or  most,  of  the  required 
solute  parameters  and  in  Table  4  are  the  b-values  used  in  equation  (9).  The 
parameters  that  were  chosen  to  be  used  in  regressions  are  as  follows: 

* 

^2  These  values  are  those  that  have  been  extensively  used  by  Kamlet 

and  co-workers ,^>2,7-10 

^2  This  is  a  trivially-calculated  parameter,  taken  as  1.0  for 

aromatic  solutes,  0.5  for  polyhalogenated  solutes,  and  zero  for 
all  others. 

These  values  were  taken  from  recent  papers  of  Kamlet  and  co- 
workers  . 

This  is  not  the  case  if  only  kovat’s  retention  indices  are  known.  It  is 
therefore  fortunate  that  Laffort  and  co-workers^  had  the  foresight  to  record 
values  of  b  for  each  of  the  five  stationary  phases  studied.  Note  that 

slightly  different  b-values  are  given  in  Laffort's  table  III.  We  are 

indebted  to  Professor  Laffort  for  suggesting  that  those  in  the  first  rwo 
would  be  the  most  suitable  to  use  in  our  equation  (9). 


0 

2^2 


log 
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f(n2) 


MR 


MR 


X 


This  is  a  new  hydrogen-bond  acidity  parameter  recently  developed 
by  Abraham  and  co-workers  using  log  K  values  for  hydrogen-bond 
coraplexation, 

These  were,  again,  taken  from  papers  of  Kamlet  and  co¬ 
workers  , 

This  Ls  a  new  hydrogen-bond  basicity  parameter,  obtained  by 
the  same  procedure  as  oc^. 

In  all  our  calculations  we  used  the  trivially  calculatable 

13 

characteristic  volume,  as  detailed  by  Abraham  and  McGowan. 

Many  of  values  were  taken  from  the  experimental  paper  of 
Abraham,  Greilier,  and  McGill.^  Other  values  have  been 

obtained  in  this  work,  see  experimental  section  and  the  section 
on  log  values. 

14  15 

Dipole  moments  were  taken  from  standard  literature  sources.  ’ 

2  2  2 

This  refractive  index  function,  defined  as  f(n  )  =  (n  -l)/(n  +2) 
was  calculated  from  literature  values  of  the  refractive  index  of 
the  solute  liquid  at  293  K  and  the  sodium-D  line. 

The  molar  refraction  was  obtained  from  the  usual  definition,  MR 
=  f(n'^).M/P,  where  M  and  P  are  the  solute  molecular  weight  and 
density;  MR  is  the  same  as  the  electron  polarisation,  P^. 

To  eliminate  the  need  for  a  new  parameter  (the  density),  a 
modified  molar  refraction  was  calculated  as  MR^  =  f(n 


In  the  present  report,  we  set  out  equations,  based  on  the  and 

values  of  Kamlet  and  co-workers ,  and  in  a  subsequent  report  we  shall 


inves  t igate 

the 

use 

of 

the  new 

parameters  cl  ^  and  1  ■  ^  . 

To  som<^  extent 

equations  (3) 

and 

(5) 

a  re 

straightforward  in  that,  apart 

from  a  and  taken 

m  m 

as  standard, 

t  he 

ot  h 

e  r 

pa  rame  t  e  rs 

anil  log  L^'’  am 

we  1  1  - de f  i  lied  sol  tit  e 

parameters. 

flowi 

f'  vr*  r 

* 

i  s  pa  r  t  1  y 

•* 

d(*rLved  from  the  soiv(M)t  paramff<*r  and 

partly  obtained  through  a  dipole  moment  correlation.  One  aim  of  the  present 

★ 

work  was  therefore  to  see  if  could  be  replaced  either  by  the  dipole  moment 
(u  or  u  )  or  some  combination  of  dipole  moment  with  f(n  )  or  MR.  Of  course, 
another  aim  is  to  obtain  the  best  regression  equation  that  could  be  used  to 
predict  new  values  of  the  chromatographi c  parameter. 

Results  and  discussion 

Regressions  were  run  for  the  relative  values  of  log  L(x)  on  the  five 
phases  listed  in  Table  1.  The  coefficients  of  the  parameters  are  listed  in  a 
series  of  Tables,  together  with  the  number  of  data  points  (n),  the  multiple 
correlation  constant  (r),  and  the  standard  deviation  (s.d.).  Also  given  are 
the  per  cent  confidence  levels  for  tlie  coefficients  (correlation  *).  The 
Tables  are  as  follows: 
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Tjble  5: 

a  ,  B  ,  Vx 

A1 1  phases 

m  m 

0 

o 

f (n")  , 

U  y  Cl  ,  B  >  V  X 

m  m 

Table  6:  ^ ,  a  ,  B  ,  1  og  ^  All  phases 

2mm 

f(n2),  a  ,  P  ,  log 

m  m 


■'),  a  ,  (3  ,  Vx 
m  m 

,  u",  a  ,  p  lo^ 
m  m 


Table  7; 


Po I y ph  ether  only 


Since  the  various  regressions  in  p,  were  not  very  good,  an  additional 
series  was  run  that  included  only  solutes  with  a  single  dominant  dipole  moment 
i.e.  solutes  such  as  acids,  esters,  nitro  compounds  and  so  forth  were 
removed.  Further  regressions  were  carried  out  as  follows: 


Table  8; 

d  2\  2 

f (n  ) ,  U  ,  a  , 

e  , 

Vx 

TCEP 

m 

ra 

f(nO,  ,  a  , 

m 

m 

log 

MR^,  a  ,  3  ,  VX 

X  mm 

MR  ,  a  ,  3  ,  log 

^  mm 

2  .  . 

Although  p  is  theoretically  a  better  parameter  to  use  Chan  p  itself,  two 
sets  of  regressions  were  run  using  p  as  a  dipolar  parameter.  In  each  set, 
regressions  were  carried  out  for  all  solutes  for  which  parameters  were 
available,  and  for  solutes  with  <5  =  0: 


Table  9: 


f(n^),  p,  am,  3m,  log 


All  phases 


Table  10: 


6 


p,  am,  3®, 
p,  am,  3m, 


log 

log 


All  phases 


A  direct  comparison  of  Tables  5  and  6  shows  quite  clearly  that 
regressions  in  log  are  always  markedly  superior  to  those  in  Vx,  and  hence 
we  shall  consider  only  the  log  correlations  henceforth.  Results  in  Table 
6  for  the  general  equation  (5)  are  quite  good,  with  correlation  constants 
ranging  from  0.9943  to  0.9805  for  the  five  phases  with  about  90  solutes.  The 
constants  in  equation  (5)  make  general  chemical  sense;  all  five  phases  are 
hyd rogen-bond  bases  and  have  no  hydrogen-bond  acidity,  and  hence  the 
coefficient  in  the  term  should  be  statistically  not  significant.  This  is 
true  for  all  the  phases  except  Znnyl  E7.  This  phase  is  supposed  to  be  a 


fluorinated  ester  of  "pyrorael  litic  acid  and  a  trihydrof  luoro  alcohol**.  The 
former  is  1 , 2 , 4 , 5- benzene  tet ra-carboxy I ic  acid,  and  it  is  possible  that  the 
commercial  product  contains  either  unesterified  carboxylic  acid  or  hydroxyl 
groups . 

One  difficulty  over  a  physicochemical  interpretation  of  equation  (5)  is 

★ 

that  polarisability  effects  are  contained  in  the  s.n^  term  as  well  as  on  the 

d.5  term.  It  would  be  particularly  useful  if  the  dipolarity  and 

★ 

pola  risabi  1  i  t  y  effects,  both  contained  in  could  be  separately  counted. 

★ 

To  this  end,  we  have  investigated  the  effect  of  replacing  the  terms  by 

2 

various  combinations  of  the  dipole  movement,  as  \x  ,  and  polarisability 

functions  such  as  f(n)  or  MR.  In  Table  6  are  results  of  a  direct  replacement 

*  2  2 

in  equation  (5)  of  and  6  by  f(n  )  and  u  .  The  overall  correlation 

constants  are  not  as  good  as  those  in  the  original  equation  (5),  but  are  not 

too  bad.  However,  most  surprisingly,  the  b.3  terras  are  statiscally 

m 

significant  for  ail  five  phases,  thus  making  the  entire  ^egression  equations 
rather  suspect  f rora  a  cheraica  1  point  of  view. 

*  .  2 

In  Table  7  are  results  of  replacing  and  6  by  either  4  and  6  or  by 

2 

u  and  MR^,  for  the  polyphenyl  ether  phase  as  an  example.  There  is  an 

excellent  correlation  with  6,  a  ,  3  ,and  log  with  r  =  0.9922  and  sd  = 

m  ra 

0.093,  which  must  be  close  to  an  exhaustive  fit.  But  once  again,  the  b.3 

m 

.  .  2  2  2 
term  is  highly  significant.  Regressions  with  f(n  )  and  u  or  MR^^  and  4  for 

TONE  using  a  restricted  set  of  solutes  that  contaifi  either  no  dipole  or  else  a 

single  dominant  dipole,  are  in  Table  8.  The  only  chemically  reasonable 

★ 

regressions  are  those  c'  the  original  form  in 

Our  conclusion  as  a  result  of  the  regressions  set  out  in  Tables  5-8  is 
that  replacement  of  the  ^^^2  symbolism  by  u"”  in  combination  with  a 
polarisability  term  leads  (i)  to  regressions  that  are  not  as  good,  and  (ii)  to 
regressions  that  contain  an  unacceptable  b.H  term. 


In  terras  of  chemical  theory,  corre  l«it  ions  of  an  ene  rgy- re  lat  ed  quantity 

such  as  log  Vq  or  log  L  with  dipole  moment  should  certainly  involve 

2 

[X  and  not  p. .  Howeve r,  we  thought  it  useful  on  an  empirical  level  to 

investigate  the  use  of  U  as  a  solvent  parameter.  Table  9  gives  details  of 

regressions  where  f(n  )  and  u  replace  6  and  n  .  Once  again,  the  b.3  term  is 

2  m 

highly  significant  for  all  five  phases.  Finally,  in  Table  10,  are  results  of 

* 

simply  making  a  direct  replacement  of  by  U  in  equation(5),  to  give: 

SP  =  C  +  s.n  +  d.5^  +  a.a^  +  ^.3^  +  l.log  (10) 

The  regressions  are  all  very  good,  and  the  only  difficulty  is  that  the 

b.3^  term  is  still  too  significant  for  the  phases  TCEP,  Polyphenyl  ether,  and 

DECS.  We  carried  out  another  set  of  regressions  using  equation  (10)  for  a 

restricted  set  of  solutes  for  which  6  =  0  (about  60-65  such  solutes).  The 

regression  coeficients  and  sd  values  are  the  best  we  have  obtained  (compare 

Table  10  with  Table  6),  and  now  the  b.B  terms  are  statistically  not 

m 

significant,  as  required  (except  for  zonyl  £7  !)-  We  are  much  encouraged  by 

the  results  in  Table  10,  and  intend  to  pursue  this  line  of  regression 

★ 

analyses.  If,  indeed,  can  be  replaced  by  Ia,  not  only  would  interpretation 

bo  much  easier,  but  it  would  be  possible  to  predict  the  dipole  parameter 
rather  easier  from  a  knowledge  of  molecular  structure. 

Work  in  Progress 

Further  work  is  in  hand  on  the  refinement  of  equation  (10)  for  the 
correlation  and  prediction  of  gas-liquid  and  gas-polymer  partition 
coe  f  f  i  c  i  en  t  s .  We  hope  that  on  the  next  report  we  shall  be  able  to  set  ovit  a 
modified  equation  (10)  that  will  deal  with  the  La f fort  data  set. 

We  4  Iso  have  in  hand  the  analysis  of  a  large  number  of  log  L  valties  (or 
rhf'  »’qui  valent  log  values)  on  the  nonpol  ir  phases  Apiezon  and  sq  ialane.  It 


will  be  possible  to  extract  from  this  data  a  rather  extended  list  of  secondary 
log  values  that  will  considerably  extend  our  data  base. 

Our  next  projected  experimentation  will  be  the  acquisition  of  log  L 

values  for  a  range  of  solutes  on  some  simple  organic  solvents,  so  that  a 
direct  comparison  can  be  made  with  GLC  stationary  phases  and  with  polymer 
phases . 

Experimental 

In  order  to  increase  the  number  of  Laffort  solutes  for  which  we  had  all 
the  parameters,  additional  log  values  were  determined,  at  the  standard 

temperature  of  298.15  K.  A  short  column  was  used  of  length  50  cm  and  internal 

diameter  2  mm,  containing  8.3*^^  w/w  of  n-hexadecane  on  chromosorb  B,  mesh  size 

45/60.  The  standard  used  was  n-octane  of  log  value  3-  677  ,  and  Log 

values  for  other  solutes  were  obtained  relative  to  n-octane,  using  a  flame 
ionisation  detector,  as  described  before."^  Results  are  given  in  Table  11. 
Attempts  were  also  made  to  obtain  log  values  for  proprionic  acid  and 

higher  carboxylic  acids,  but  without  success,  but  further  attempts  will  be 
made  to  determine  these  quantities,  either  directly  or  indirectly. 
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Table  A 


J  in  u.i  t  i  on  (  M  )  a  nd  (  V  )  . 


Vdluns  ol  b  nso 

Stationary  phase  b 

Zony I  0.203 

CWAX  1540  0.214 

TCEP  0.176 

Polyph  ether  0.2b2 
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Table  11* 


Further  values  of  log  obtained  at  298.15  K 

relative  to  n-octane  standard. 


Solute 

log  L 

n-Oc  tane 

3.677 

Benzoni t  ri le 

3.994 

Benza Idehyde 

3.985 

2 -Methyl -3-pentanol 

3.  183 

3 -Methyl -3-pentanol 

3.227 

3“Hexanol 

3.440 

Ethyleneglycol  diacetate 

4.083 

Pyrrole 

2.866 

2-Ch lorophenol 

4.937 

N-Met  hy 1 ani 1 ine 

4.494 

2 -Ch 1 orot o 1 uene 

4.  160 

Methyl  benzoate 

4.634 

Mes i t y 1 ene 


4.399 


